Probing nearest-neighbor correlations of ultracold fermions in an optical lattice 
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We demonstrate a probe for nearest-neighbor correlations of fermionic quantum gases in optical 
lattices. It gives access to spin and density configurations of adjacent sites and relies on creating 
additional doubly occupied sites by perturbative lattice modulation. The measured correlations 
for different lattice temperatures are in good agreement with an ab initio calculation without any 
fitting parameters. This probe opens new prospects for studying the approach to magnetically 
ordered phases. 
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The Fermi-Hubbard Hamiltonian is one of the central 
models for understanding strongly correlated electron 
systems in condensed matter physics. It incorporates in- 
triguing phenomena such as Mott-insulating behavior or 
spin ordered phases and is a prominent candidate for de- 
scribing the origin of high Tc superconductivity. Despite 
tremendous theoretical effort, several questions still re- 
main open, in particular concerning the low temperature 
phases. Here, ultracold atomic gases trapped in optical 
lattices offer the opportunity to address these questions 
in a very clean way, as they constitute an almost ideal im- 
plementation of the Hubbard model jj 2n ■ The recent re- 
alization of a fermionic Mott insulator [5J|3] demonstrates 
the unique tunability and control of these systems. 

Considerable experimental efforts are currently di- 
rected towards reaching the low-temperature regime of 
quantum magnetism in a two-component quantum gas. 
Detection of the antiferromagnetically ordered state has 
been proposed via noise correlation or Bragg scattering 
measurements IS] • These observables reveal long range 
spin ordering and consequently only show a strong sig- 
nature well below the critical temperature. However, a 
probe sensitive to local magnetic correlations [3 [S] for 
studying the approach to magnetic ordering of fermions 
close to the transition point is missing so far. 

In this Letter we demonstrate a simple method for 
probing the nearest-neighbor correlations of strongly in- 
teracting repulsive fermionic gases in optical lattices. 
The correlation function is given by 

'Pi.i+i = (1 - Ui,^) (1 - Ui+i^^)), (1) 

where a = {f, 4-}, are opposite spins and i, adjacent 
sites. This probe determines the probability of finding 
singly occupied neighboring sites with opposite spins. 

The experimental strategy for detecting the correlator 
Vi^i+i relies on exciting the system by a periodic modula- 
tion of the lattice depth [5]. The corresponding modula- 
tion in kinetic energy leads to tunneling of particles to ad- 
jacent sites. If two particles of opposite spin are located 
on neighboring sites, additional double occupancies (dou- 
blons) are created as shown in Fig. IT] The resulting dou- 



blon production rate is sensitive on the nearest-neighbor 
density and spin correlator Vi^i+i. In the following we 
describe and characterize this experimental method and 
show that in the perturbative regime the frequency inte- 
grated doublon production rate is given by Vi^i+i [lOHIl]. 
We then use this technique to measure nearest-neighbor 
correlations as a function of temperature covering the 
regime from a paramagnetic Mott insulator to a strongly 
interacting metallic state. The results are in good agree- 
ment with the predictions of an ah initio theory without 
any fitting parameters. 

The experimental sequence used to produce a quan- 
tum degenerate Fermi gas has been described in detail 
in previous work |14L I15j . In brief, a balanced spin mix- 
ture of ''''K atoms in the mp = —9/2 and —5/2 magnetic 
sublevels of the F = 9/2 hyperfine manifold is evapora- 
tively cooled in an optical dipole trap. For samples of 
N = 80(7) X 10"^ atoms we reach temperatures as low as 
14% of the Fermi temperature Tp. Subsequently we ramp 
up a three-dimensional optical lattice of simple cubic ge- 
ometry and lattice constant d = 532 nm. The lattice 
depth is increased in 0.2 s to final values of 7.0(7) or 
10(1) £'r , where i?R — h'^ /8mcP is the recoil energy, h is 
Planck's constant and m denotes the mass of '^^K. The 
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FIG. 1: (color online) Probing nearest-neighbor correlations 
for different phases. A periodic lattice modulation causes tun- 
neling of particles to neighboring sites. The number of created 
doublons strongly depends on the state of the many-body sys- 
tem (strongly interacting metal, paramagnetic Mott insulator 
or antiferromagnet) and can be used to determine the nearest- 
neighbor correlator Pi,i+i [TO] . 
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hopping t is inferred from Wannier functions [T5] and 
the on-site interaction energy U is obtained from lattice 
modulation spectroscopy [3]. The underlying trapping 
potential has a mean frequency of w/27r = 70.1(5) Hz for 
7 En and w/27r = 77.3(7) Hz for 10£:r. With this pro- 
cedure we create samples where the core is in the Mott 
insulating regime [15j . 

After this preparation, the lattice depth is modu- 
lated by SV along all three axes in time t according 
to V{t) — V + (51^sin(27rj/r). This results in a mod- 
ulation of both the hopping and the on-site interaction 
with amplitudes 6t and SU respectively, creating addi- 
tional doubly occupied sites as compared to the initial 
state. The increase in the number of doublons is maxi- 
mal when the modulation frequency coincides with the 
doublon energy (resonant excitation at v — U/h). Af- 
ter the lattice modulation, the fraction of atoms on dou- 
bly occupied sites D = '^'^iini^^rii^i) /N is measured by 
mapping doublons into a different spin state, subsequent 
Stern-Gerlach separation and absorption imaging [3 . 

Fig. [2] (a) shows the evolution of double occupancy as 
a function of lattice modulation time for resonant excita- 
tion. After a steep initial rise, we observe a saturation of 
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FIG. 2: Evolution of double occupancy as a function of the 
lattice modulation time r for resonant excitation. The lattice 
depth is set to 7 E-r, {U/6t — 4.1) and the modulation strength 
is 5V/V = 0.1. (a) The induced double occupancy saturates 
for large times, which is well captured by an exponential fit 
(solid line), (b) At low modulation times D increases linearly, 
from which the doublon production rate F is obtained by a 
linear fit. The lattice depth was set to 10 -Br for this measure- 
ment {U/6t = 10.6). (c) On the timescale of the modulation 
period hjU the double occupancy shows an underlying low 
amplitude sinusoidal modulation. The solid line is a fit with 
fixed frequency U jh. For clarity, &V jV was increased to 0.2. 
Error bars in double occupancy denote statistical errors from 
multiple measurements. 



the induced double occupancy on a timescale on the order 
of the tunneling time hjt. The saturation value depends 
only weakly on the modulation strength and reaches typ- 
ical values of 20 — 30%. In contrast to previous work, 
where only this saturation regime was considered [31 [13] , 
our high accuracy in the determination of double occu- 
pancy allows us to perform measurements in the weak 
excitation limit. Here we find that double occupancy 
increases linearly with time as shown in Fig. [2j[b) . We 
extract the normalized doublon production rate F from 
the slope AD/ At of a linear fit to the data 
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On shorter timescales an underlying oscillatory response 
at the modulation frequency v is observed, Fig.[2](c). 

The experiment can be well understood in the frame- 
work of time-dependent perturbation theory, which we 
outline below. The main result is that the frequency 
integrated doublon production rate is proportional to 
the nearest-neighbor correlator T'i^i+i. The Hamiltonian 
of the system can be written as H — Hq + ffpcrt(''")j 
where i?o = —tHt + UHu is the Fermi-Hubbard Hamil- 
tonian and the time-dependent perturbation is given by 
^^port(''') = {5tHt + 5UH\j)sui{2ttvt). Assuming that 
the doublon production rate is equivalent to the total 
energy absorption rate (which becomes exact at half- 
filling and [/ 3> t), the perturbation in U can be mapped 
to an increased tunneling perturbation with amplitude 
5t/t = St/t — SU/U [17 . The response of the system to 
second order in perturbation theory is then given by 

Dir) = D(0) + [ y j x^'Hi') sin(2^z.T) 



+ 



X^^^ (v) + osc. terms 



(3) 



in good agreement with the experimental observations of 
Fig. [21(b), (c) using F = x^^\i')iSt/t)^. The susceptibili- 
ties x^^n^) x'^'^H^) correspond to the linear response 
and the non-oscillatory part of the quadratic response re- 
spectively. This is recovered by Fermi's Golden Rule 
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Here |0) and \n) denote the unperturbed and excited 
states of _ffo, is their energy difference and bHu 

counts the number of additionally created doublons [11] . 
Evaluating all contributing matrix elements in the limit 
U ^ t, the sum over all excited states is then equal to 
Vi^i+i- To compare with experiments we consider the 
normalized frequency integrated response 



i? = ^ y diyT{v) = 27r2 




(4) 
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where V = ^i'Pi^i+i/N is the system averaged corre- 
lator and z the connectivity of the lattice. In the per- 
turbative regime R therefore gives direct access to the 
nearest-neighbor correlator V of the unperturbed initial 
state. 

We validate that the experiments are performed in the 
weak excitation regime by studying the scaling of R with 
the relative modulation amplitude 6t/t [2^. We measure 
the doublon production rate as a function of the modu- 
lation frequency, from which the modulation spectra in 
Fig. [3] are obtained. The frequency integrated response 
is then determined by a gaussian fit to each spectrum. 
The result is plotted as a function of the lattice modula- 
tion amplitude St/t (inset of Fig. We find a scaling 
exponent of 2.09(18), in very good agreement with the 
expected value of 2 predicted by second order pertur- 
bation theory. We can thus infer the nearest-neighbor 
correlator V from the frequency integrated response. 

Further information can be obtained from the line- 
shape of the modulation spectra, which reveals the den- 
sity of states of the excitations. At half filling and tem- 
peratures well above the Neel transition, the density of 
states has an approximately triangular shape of full width 
~ 3zt [12]. However, the trapping potential is expected 
to broaden the spectrum and introduce deviations to the 
lineshape [lOj . This is consistent with the experimental 
data, which is well captured by gaussian fits with 1/e^ 
diameters of Azt. The doublon production rate on reso- 
nance shows the same scaling behavior as the integrated 
response R, with a scaling exponent of 1.98(11). This 
allows us to determine V from the resonant doublon pro- 
duction rate alone assuming a gaussian density of states 
as in Fig. |3]pg. 

We now use doublon production rate measurements to 
determine the nearest-neighbor correlator 7^ as a function 
of entropy for U/6t = 4.1. We prepare samples with dif- 
ferent entropies per particle by adding a variable waiting 
time in the optical dipole trap of up to 2 s, which results 
in heating due to inelastic scattering processes. The en- 
tropy per particle Sin before loading into the lattice is 
inferred from Fermi fits to the momentum distribution of 
the cloud after expansion. This is a lower bound for the 
specific entropy in the lattice, as non-adiabatic processes 
take place during the loading. An upper bound Sout is 
given by the entropy measured after reversing the loading 
procedure. 

The nearest-neighbor correlator rapidly decreases with 
increasing entropy, as shown in Fig. |4](a). This behavior 
has a simple physical interpretation: in a harmonically 
trapped lattice system higher temperatures lead to an 
increased cloud size, which results in a large number of 
empty sites. The probability of finding two neighboring 
singly occupied sites is therefore strongly reduced. This 
qualitative picture is confirmed by ab initio calculations 
of the nearest-neighbor correlator. We use a high tem- 
perature series expansion (HTSE) up to second order in 
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FIG. 3: (color online) Doublon production rate F as a function 
of the lattice modulation frequency v, measured for different 
modulation amplitudes St/t. The experiments are performed 
at U /6t = 10.6 and V = IO-Er. The shaded areas are gaussian 
fits to the spectra, which are used to extract the frequency 
integrated response R. The two vertical dashed lines denote 
twice the three-dimensional bandwidth Azt (with z = 6). The 
inset is a double-logarithmic plot of R for various modulation 
amplitudes, where the dashed line is a power law fit. Error 
bars denote the fit errors. 

t/k^T, where fee is the Boltzmann constant [THllin]- Due 
to the harmonic confinement, the trap averaged correla- 
tor V needs to be evaluated. This is done using a local 
density approximation, which is an excellent assumption 
in this temperature regime |20j . The system parameters 
are calibrated by independent methods I15j. The results 
of this calculation are shown in Fig. |4|^a). 

We find quantitative agreement between the measured 
nearest-neighbor correlator and the theoretical predic- 
tions without any fitting parameters. The theoretical 
entropies s corresponding to the measured values of V 
lie in between the lower and upper experimental bounds 
Sin and Sout • Nearest-neighbor correlation measurements 
thus allow us to determine entropy and temperature in 
the lattice for regimes where theory is still reliable and a 
quantitative comparison is possible. In contrast to ther- 
mometry in the optical dipole trap, the correlator V is 
a direct observable in the lattice and does not rely on 
adiabaticity assumptions during the loading procedure. 
A comparison of these two methods suggests increased 
heating during the loading of the lattice for colder initial 
temperatures. 

Further insight can be obtained from the theoretical 
model by investigating the spatial distribution of the 
nearest-neighbor correlator over the trap. From the in- 
ferred entropies s we calculate the profiles shown in Fig. 
|4|^b) and (c) for the parameters of our system [26]. At 
half filling and deep in the Mott insulating regime Vi,i+i 
is expected to be close to 0.5, whereas thermal excita- 
tions reduce this number in the metallic phase. This is 
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FIG. 4: (color online) Measurement of the nearest-neighbor 
correlator V. (a) Dependence on entropy per particle in a 
trapped lattice system for U/6t = 4.1, 1/ = 7 En and 5t/t = 
0.26. Solid and open circles denote lower and upper bounds 
Sin and Sout for the entropy per particle in the lattice for each 
measured value of V. The black solid line is the calculated 
nearest-neighbor correlator obtained from second order HTSE 
without any fitting parameters. The entropy s inferred from 
comparing each measured value of V to theory lies in between 
the experimental bounds. For clarity the correspoding lattice 
temperatures T and measured doublon production rates F are 
also included. Error bars denote statistical errors from several 
measurements, (b) and (c) show the calculated distribution 
of Vi,i+i in the trap for the lowest and highest entropies Scom 
and Shot , while the area under these curves corresponds to the 
measured correlator V. 



confirmed by the values in the core of the system, with 
0.44 for the coldest and 0.17 for the hottest point. The 
decrease in Vi^i+i thus signals the transition from a para- 
magnetic Mott insulator to a strongly interacting metal. 
As the onset of local spin correlations corresponds to an 
increase of Vi^i+i above 0.5, nearest-neighbor correlation 
measurements are a promising tool for studying the ap- 
proach to the antiferromagnetic phase [22] . 

In conclusion, we have measured nearest-neighbor cor- 
relations of ultracold fermions in optical lattices by de- 
termining the response of the system to a weak lattice 
modulation. This observable is well suited for thermom- 
etry in the lattice and can be used to explore novel cool- 
ing schemes |21j . The technique opens new prospects for 
studying the approach to the antiferromagnetic phase, 
since the regime between a paramagnetic Mott insula- 
tor and an antiferromagnet is governed by the forma- 
tion of short-range magnetic correlations. In the future. 



nearest-neighbor correlation measurements might give in- 
sight into resonating valence bond ground states, where 
singlet correlations on neighboring sites are expected to 
occur in the absence of long range ordering [23] . 
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